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Abstract 

We present a calculation of the generalized parton distributions (GPDs) of the photon when the 
helicity of the initial photon is different from the final photon. We calculate the GPDs using overlaps of 
photon light-front wave functions (LFWFs) at leading order in electromagnetic coupling a and zeroth 
order in the strong coupling a s , when the momentum transfer is purely in the transverse direction. 
These involve a contribution of orbital angular momentum of two units in the LFWFs. We express 
these GPDs in the impact parameter space. 
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Introduction 



Generated parton distributions (GPDs) of the nocleon are unified objects giving a w ld e 
range of information on nuclear structure and spin [lj. These are non-perturbative objects 
appearing in the factorized amplitude of exclusive processes like deeply virtual Compton scat- 
tering (DVCS) and meson production; and can be expressed as an off-forward matrix element 
of light-cone bilocal operators. In [2| the amplitude of the DVCS process on a photon target 
l*{Q)l ~^ 77 a t high Q 2 is written in terms of photon GPDs. These photon GPDs were cal- 
culated at leading order in electromagnetic coupling a and zeroth order in the strong coupling 
a s and upto leading logs; in the kinematical limit that there is no momentum transfer in the 
transverse direction. In fact the parton content of the photon is known to play an important 
role in high energy scattering processes. The parton distributions of the photon are now well 
understood both theoretically and experimentally [3j. On the other hand, the GPDs and gen- 
eralized distribution amplitudes (GDAs) of the photon J4] are much less investigated objects. 
In a couple of recent works 0, Q], we extended the calculation of photon GPDs in the more 
general kinematics when the momentum transfer has both transverse and longitudinal compo- 
nents. We have developed an overlap representation using the light-front wave function of the 
photon. We also showed that the impact parameter space interpretation of the photon GPDs 
give a 3D position space description of them. In another recent work [7j], GPDs of the photon 
have been used to investigate analyticity properties of DVCS amplitudes and related sum rules 
for the GPDs. 

As we know, in the DVCS process eP — > e'yP, the helicity of the proton may or may 
not flip due to the scattering. When the proton helicity is flipped, the DVCS amplitude is 
parametrized in terms of the GPD E [lj]. This flip requires non-zero orbital angular momentum 
in the overlapping light-front wave functions (LFWFs) and is not possible unless there is non- 
zero momentum transfer in the transverse direction. For a transversely polarized nucleon, this 
gives a distortion of the parton distributions in the transverse position or impact parameter 
space jsj]. In two previous articles, we calculated the impact parameter space representations 
of the photon GPDs when the helicity of the photon is not flipped. In this work, we calculate 
the GPDs that involve helicity flip of the photon and represent them in impact parameter 
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space. Like the proton, these involve overlaps of LFWFs of the photon, with non-zero orbital 
angular momentum (OAM). The corresponding parton distributions in the impact parameter 
space show distortions related to the orbital angular momentum of the LFWFs. 



GPDs of the photon with helicity flip 



The GPDs of the photon can be expressed as the following off-forward matrix elements 

BB 



here | 7(-P), A) is the (real) photon target state of momentum P and helicity A. We work in the 
light-front gauge A + = 0. We use the standard LF coordinates P^ 1 = P° ± P 3 , y^- = y° ± y 3 . 
Since the target photon is on-shell, P + P~ — P ±2 = 0, the momenta of the initial and final 
photon in the most general case of momentum transfer are given by: 

P = ( P+ , 1 , ) , (3) 

The four-momentum transfer from the target is 

A = P-P' = ((P + ,A\ tJ ^Pj , (5) 

where t = A 2 and £ is called the skewness variable. In addition, overall energy-momentum 
conservation requires A~ = P~ — P'~, which connects A ±2 , and t according to 

(1 - C)t = -A ±2 . (6) 

In order to calculate the above matrix element, we use the Fock space expansion of the 
photon state, which can be written as Q 

| 7 (P), A) = VN[<f(P, A) | 0) + j i dk ^ J {dh}^2(2n) 3 P+6 3 (P - h - k 2 ) 



fa(h, k 2 , cti, a 2 )b\h, a 1 )d\k 2 , a 2 ) | 0) (7) 

where \/N is the normalization of the state; which in our calculation we can take as unity 
as any correction to it contributes at higher order in a. {dk} = J ^f=p=, 02 is the two- 
particle (qq) light-front wave function (LFWF) and <j\ and o 2 are the helicities of the quark 

k + 

and antiquark. The wave function can be expressed in terms of Jacobi momenta Xi = -p+ and 
qf- = ki~ — XiP^. These obey the relations Yli x i = 1, Tlill = 0- The Lorentz boost invariant 
two-particle LFWFs are given by ijj 2 (xi,q^) = 4> 2 \J P + . vjj 2 (xi,q^) can be calculated order by 
order in perturbation theory. The two-particle LFWFs for the photon are given by 



W2si,S2\ X 1 ( l > = 2~t ±\2 / Xsi 

2 m i +(q^) 3 



x(l—x) 



2(2tt) 



X 



\ \p _g . m ^ 

-a — i— ra 

1 — x x(l — X) 



X-sA* (8) 



where m is the mass of q(q). A is the helicity of the photon and si,s 2 are the helicities of 
the q and q respectively. We have used the two-component form of light-cone field theory 9] , 
namely the component A~ of the photon field is constrained in the gauge A + = and can be 
eliminated from the theory. So one has only the transverse components of the photon field A- 1 . 
Likewise, the 'bad' component of the fermion field ip^ is eliminated using constraint equation 
and ip^ is written in terms of two-component spinors, Xs 

The GPDs can be written in terms of the overlaps of the LFWFs as follows : 



F q = j d 2 q ± dx 1 5(x - Xi)^ 2 * A (x u q ± -(l- x^A^i/fifa, q L ) 

- J d 2 q ± dx 1 8(l + x - xi)il)* 2 x ' {xi, q L + XiA- 1 )^ q L ) (9) 

We calculate the photon GPDs using overlaps of light-front wave functions. We take the 
momentum transfer to be purely in the transverse direction, unlike |2|, where the momentum 
transfer was taken purely in the light-cone (plus) direction. GPDs in this case can be expressed 
in terms of diagonal (particle number conserving) overlaps of LFWFs. When there is non- 
zero momentum transfer in the longitudinal direction, there are off-diagonal particle number 



changing overlaps as well, similar to the proton GPDs 10] . 
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The transverse polarization vector of the photon can be written as : 

1 



y/2 



(=Fl,-<) 



(10) 



We extract the GPD that involves a helicity flip of the target photon from the non- vanishing 
coefficient of the combination (e^e^ + e^e 2 ^). The corresponding GPD without a helicity 
flip of the photon contains a leading logarithmic term at leading order in a and zeroth order in 

nn 

strong coupling constant and has been discussed in two previous articles |5|,|6[. The GPD with 
helicity flip is given by : 

-,2 
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x(l — x) Ji — (1 — x)Ii 



'IV, 



The integrals I\ and I 2 are given by : 
I x = [ d 2 q 



2 Atf) 2 -(<?)') 



j D X D 2 J * D X D 2 

where q l and q 2 are the x and y components of q L and A 1 and A 2 are the x and y components 
of A ± respectively. The denominators are given by : 



D\ = (g -1 ) — m 2 x(l — x) + m 2 
D 2 = (q^f + (1 - x) 2 (A x ) 2 - 2q ± ■ A ± (l - x) - m 2 x(l - x) + m 2 . 



In order to simplify the above expression we use the formula 11 1 

1 1 ro ° 



a k - l e~ aA da. 



A k T(k) Jo 
The integrals can be written in the form : 



(12) 



h = ((A 1 ) 2 - (A 2 )V(1 - x) 2 J 1 dq J 2 = ((A i) 2 _ (A 2 ) 2 )tt(1 - x) J 1 dq A 



(13) 



9). 



where 



B(q) = m 2 (l- x{l - xj\ + q(l - q){\ - x) 2 (A L f. 



(14) 



So we have; 
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(a) 
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x = 0.1 
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FIG. 1: (Color online) Plots of E\{x, A ) vs A , A for fixed values of x and different t. t is in GeV 



The innermost surface is for the smallest value of —t. 
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1 dq 
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((l-9) 2 -(l-g)) 



(15) 



The above has the expected quadrupole structure coming from ((A 1 )' - (A 2 ) ). As the pho- 
ton is a spin one particle, in order to flip its helicity, the overlapping light-front wave functions 
should have a difference of orbital angular momentum of two units, which manifests itself in 
the quadrupole structure. This is in accordance with a similar observation for the helicity-flip 
GPD E for the proton, which needs overlapping LFWFs of orbital angular momentum ±1 unit 
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From the off-forward matrix element F q we extract the GPDs that flip the helicity of the 
photon by calculating the coefficient of the combination (e^e 2 ^ + e^e^) which gives: 



4tt 2 



(x - x') h - (1 - x'% 



h = 7T I dq ( (Af - A 



21 B(q) 
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FIG. 2: (Color online) Plots of (a) Ei(x,A- L ) vs x for different values of t in GeV 2 and (b) qi(x,b' i 

and /3 = tan _1 p 



vs. x for different values of b in GeV 1 . <ft = tan and /3 = tan 1 
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1 dq 

W) 



-q)(l-q) 



(16) 



Here x' is the longitudinal momentum fracion of the quark in the final photon LFWF. As 
we have taken the momentum transfer to be purely in the transverse direction, x' = x and 
Ei = 0. Ei and E\ are only two independent structures that cause helicity flip of the photon. 
Other helicity-flip GPDs that can be constructed from other combinations of the polarization 
vectors can be related to these by phase change in the A 1 - plane. However a proper counting 
of the photon GPD (both helicty non-flip and helicity flip) can only be done in a formal 
parametrization of Eqs. (JT|) and (f2]). 



Like the GPD E of a spin 1/2 particle for example a dressed electron/quark 13|, the helicity 
flip photon GPD has no logarithmic term depending on the hard scale of the process Q 2 , which 
is the virtuality of the probing photon. Starting from the expressions of photon GPDs, we 
define the parton distributions s| with the helicity flip of the photon in transverse impact 
parameter space as: 
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(c)l 1 (d)L 

FIG. 3: (Color online) Plots of qi(x, b- 1 ) vs. b 1 , b 2 for different values of A max . b 1 and b 2 are in GeV -1 
and A. max is in GeV. x = 0.3. 



qi {x^) = ^- 2 j d 2 A ± e- iAX ' b± E 1 (x,A ± ); (17) 
where t = — (A -1 ) 2 and b 1 - is the transverse impact parameter conjugate to A -1 . One then gets 
q 1 (x,b ± ) = -LJ d 2 A^e~ ih± ^{{A l f - (A 2 ) 2 )/(x)Q(x,t), (18) 

where 

ae 2 C 1 rln 

f(x) = - x) 3 , Q(x, t) = J ^ ((i - qf - (i - *)); (19) 

This can be written as, 
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qi{x,b L ) 



d 2 



d 



An 2 \d(b 2 Y d(b l 
1 / d 2 d 2 



^) J d 2 A^e-* b± - A± f(x)Q(x,t) 

AdA J (bA) f(x)Q(x,t) 



(20) 



2n\d(b 2 ) 2 d{b l f 

Here A =| A L | and b —\ b |. Using the integral representation for the Bessel function Jq(x 
the above can be written as, 



qi(x,b 



We then get 



1 / d 2 



d 2 



2n\d(b 2 Y d(b l Y 



AdA 



i r 



TT 



cos(6A - sm8)d6 f(x)Q(x,t) . (21) 



i poo i rn 

q 1 (x,b ± ) = — AdA- (P 2 (b,A,9)-P 1 (b,A,e))d9 f(x)Q(x,t); 



where 



2tt 



TT 



(22) 



P 2 (b,A,6) = -^Asin^ (6 2 ) 2 6Acos^Asin^sin^ + 



(tffsmlbAsmO 



Pi(b,A,9) = -^Asin^ (fe^^Acos^Asin^sin^ + 



(6 2 ) 2 sin(6Asin^ 



(23) 



(24) 



Numerical Results 



We next discuss the numerical results. In Fig. 1, we have shown the helicity flip GPD of the 
photon as functions of A -1 . As we mentioned before, we took ( = 0. In this kinematical limit, 
the GPDs are represented by overlaps of two-particle photon LFWFs. When ( is non-zero or 
the momentum transfer between the initial to final photon has a longitudinal component, off- 
diagonal particle number changing overlaps of LFWFs has to be considered as well 10j. We take 
the mass of the quark and the antiquark in the photon to be the same and equal to 3.3 MeV. As 
seen in Eq. (Q, the GPDs have two contributions. When x is positive, the contribution comes 
from the active quark in the photon; and when x is negative, the active antiquark contributes. 
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FIG. 4: (Color online) Plots of qi(x, 6 ) vs 6,6 for different values of a; and at A maa; = 4.0 GeV. 



In the numerical plots, we take < x < 1 and show the quark contribution. Fig. 1 shows the 
helicity flip GPD E\{x, A ) as functions of A 1 and A 2 and for different values of x for fixed 
values of t. Ei(x, A^) is zero when A 1 = A 2 . The curvature is sharper as | t | decreases. As we 
already saw, Ei(x, A -1 ) has a quadrupole structure in A 1 - plane coming from the (A 1 ) -(A 2 ) . 
Such quadrupole structure is due to the spin flip of a spin one particle and corresponds to 
an overlapping LFWF with two units of orbital angular momentum. The structure can be 
contrasted with the GPD E(x, A -1 ) of a spin 1/2 composite particle like a dressed electron or 



a proton 



13Nl5|. It is to be noted that the off- forward matrix element similar to Eq. ([I]) for 



a proton target is parametrized in terms of the GPDs H and E. When the final proton has 
the same helicity as the initial proton, and ( is non-zero, both the GPDs contribute. However 
when the helicity of the proton is flipped then only the GPD E contributes. A parametrization 
of the off- forward matrix element for a spin one massive target like deuteron was given in 
So far no such parametrization is available for the photon GPDs and in this work as well as in 
two previous publications 5|, |6[ we calculate the full off-forward matrix element using overlaps 
of photon LFWFs. In Fig. 2(a) we have plotted the helicity flip photon GPD Ei(x, A -1 ) vs. x 
for different values of t and a fixed value of <p = tan -1 ^-. As before, we have plotted for the 
region < x < 1, where the contribution to the GPD comes from the active quark. The peak 
of Ei(x, A -1 ) increases as —t increases and also shifts towards larger value of x. The GPD is 
zero both at x = and x = 1. In fact, the GPD is zero when A x = 0. This is because in 
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order to flip the helicity one needs non-zero OAM in the two-particle LFWFs and the OAM 
is zero when there is no momentum transfer in the transverse direction. At x = and x = 1 
all momenta are carried by either the quark or the antiquark in the photon. Then there is no 
relative motion and no OAM contribution. In Fig. 2 (b) we have plotted the Fourier transform 
(FT) of the helicity-flip photon GPD, g 1 (x,6 _L ) vs. x for different b =\ b 1 - | at a fixed value 
of f3 = tan -1 ^. gi(x,6 _L ) is symmetric with respect to x = and x = 1 and maximum when 
x = 0.5, that is when the quark and the antiquark carry equal momenta. As seen in Eq. (|2ip . 
qi(x, b^) has a quadrupole structure, that comes because it involves a helicity flip of a spin one 
object (photon). This quadrupole structure is visible in the 3D plots of Figs 3 and 4. In the 
ideal definition of the Fourier transform, the limits of the A x integration should be from to 
oo. As we saw for the photon GPDs that do not involve a helicity flip , the A 1 - independent 
terms then give a <5(& _L ) in impact parameter space. For non-zero A -1 , we get a smearing in b^. 
For the GPD with helicity flip, from Eq. (l2Tj) . we see that it involves a distortion in b x space. 
The GPD as well as its FT is zero when A -1 = 0, which means that it is purely an effect of the 
orbital angular momentum of the LFWF. In the actual numerical calculation we have imposed 
an upper limit on the A -1 integration, denoted by A max . Fig. 3 shows a plot of qi(x, b^) vs. 
b 1 and b 2 for a fixed value of x = 0.3 and different values of A max . It is seen that as A max 
increases the peaks become sharper, which means that the distortion in b 1 - space moves closer 
to the origin. Fig. 4 shows plots of qx(x,b ) vs. b 1 and b 2 for a fixed value of A max and two 
different values of x. The magnitude of the peaks depend on x. 

Conclusion 

In this work we have calculated the GPDs of the photon when the helicity of the target 
photon is flipped. We expressed the GPDs in terms of overlaps of photon LFWFs. In the 
kinematics when the momentum transfer between the initial and the final photon is purely in 
the transverse direction, the GPDs involve diagonal overlaps of two-particle LFWFs at leading 
order in the electromagnetic coupling a and zeroth order in the strong coupling a s . Such two 
particle LFWFs of the photon can be calculated in light-front Hamiltonian perturbation theory. 
Taking a Fourier transform of the GPDs with respect to A 1 - we obtained impact parameter 
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dependent parton distributions. Like the proton GPD E the helicity flip GPD of the photon 
represents a distortion of the parton distribution in the impact parameter space. This is due 
to the orbital angular momentum contribution coming from the LFWFs. As photon is a spin 
one object, one needs OAM of two units in the overlapping LFWF to flip the helicity. The 
expected quadrupole structure is visible in impact parameter space. For the proton, such 
distortion in b 1 - space has been found to be related to the Sivers function in some models. It 
will be interesting to check if such relations exist also for the photon. For this it is necessary 
to have a parametrization of the off-forward as well as the transverse momentum dependent 
matrix elements for the photon. 
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